Focusing in Multiwell Potentials: Applications to Ion Channels 
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We investigate out of equilibrium stationary distributions induced by a stochastic dichotomous 
noise on double and multi-well models for ion channels. Ion-channel dynamics is analyzed both 
through over-damped Langevin equations and master equations. As a consequence of the external 
stochastic noise, we prove a non trivial focusing effect, namely the probability distribution is con- 
centrated only on one state of the multi-well model. We also show that this focusing effect , which 
occurs at physiological conditions, cannot be predicted by a simple master equation approach. 

PACS numbers: 05.50.+q, 75.10.Hk, 75.10.Pq 



I. INTRODUCTION 

Interaction of physical, chemical or biological systems 
with fluctuating stimuli, both random and/or periodic, 
has in recent years received a lot of attention in a vari- 
ety of contexts. It gives rise to new effects in nonlinear 
systems that cannot be observed under stationary per- 
turbation. One of the most studied examples of such 
phenomena is the stochastic resonance [THE] ■ It is an ef- 
fect when a noise with suitable properties, either intrinsic 
or added to the system, improves the system response to 
weak time-dependent signals. The signal-to-noise ratio 
can paradoxically be boosted by increasing the level of 
noise present in the system. The SR was first suggested 
as an explanation for the cyclic ice ages observed in the 
earths climate [T] but has been since expanded as a plau- 
sible explanation of various effects in other areas. Of 
particular interest seems to be the application to explain 
the detection of very weak signals in noisy environments 
observed in biological systems, e.g. ion channels [HI El E] - 

Another example is the dynamics of ratchets. The 
term refers to spatially asymmetric potentials which, in- 
teracting with Brownian fluctuations generate a directed 
drift in the system. The scale of the effect can be con- 
trolled by adjusting the properties of the noise. A simi- 
lar ratchet-like effect has been observed in spatially sym- 
metric potentials subject to temporally asymmetric noise 
(so-called correlation ratchets) [9]. 

An interesting example of such fluctuation-induced 
phenomena is called the resonant activation [TU1 E] . It 
is an effect where the diffusion over a fluctuating poten- 
tial barrier is correlated with the fluctuation rate and a 
resonant effect is observed. As noted in [TU] if the barrier 
fluctuations are slow compared to the natural time scale 
of barrier crossing then the effect can be described using 
an adiabatic approximation. For very fast oscillations 
the barrier crossing rate would be that for the average 
barrier height. For intermediate barrier fluctuation rates 
a resonancc-likc enhancement of barrier crossing rate is 
observed. This has been applied e.g. to chemical reaction 



rates [T^] and ion transport through channels [15] . 

A related effect, and the one we investigate in this 
paper is the nonequilibrium kinetic focusing first pro- 
posed in [14] . It is an effect where a particle moving in 
a ratchet-like potential subject to external dichotomous 
noise fluctuations gets trapped in one of the energy wells. 
This mechanism was applied to gating kinetics of voltage- 
gated ion channels, achieving a focusing of the channels 
into a particular conformational state corresponding to 
this energy well. First experimental study of this effect 
in Shaker K+ channels was reported in [15] . The phe- 
nomenon of particle trapping in correlation ratchet was 
also investigated in [5]. 

Ion channels are membrane proteins in biological cells 
that form gated channels for a controlled exchange of 
physiologically important ions, such as sodium or potas- 
sium, between the cytosol and the extracellular medium 
[16] . They play a very important role in various phys- 
iological processes (nerve impulses, muscle contraction 
etc.), and several human and animal disorders have been 
linked to malfunctioning ion channels. Therefore there 
have been numerous studies aimed at better understand- 
ing of channel gating kinetics and ultimately controlling 
it. Kinetic focusing of ion channels seems to be very 
promising for both of these goals. 

Gating kinetics of voltage-gated ion channels is a reflec- 
tion of conformational changes of the channel molecule. 
From a physico-chemical standpoint the gating process 
has been modeled as a particle (so-called gating parti- 
cle) moving in a certain multi-well energy landscape (the 
energy of molecular conformational states). Mathemat- 
ically, the most commonly used description is a discrete 
Markov chain, where different Markov states correspond 
to minima in the energy profile and the transition rates 
between the states reflect the height and width of the 
energy barrier. The time evolution of the system is de- 
scribed by the master equation. These types of mod- 
els dominate the literature about ion channels and have 
been reasonably successful in explaining physiologically 
relevant channel effects (see e.g. [17J. However, they are 
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coarse approximations only, completely disregard intra- 
well motions of the gating particle, and work only if the 
energy barriers separating the wells are sufficiently high. 
More accurately the process can be described by the over- 
damped Langevin equation [T]. 

In this paper we study both approaches to simple mod- 
els of ion channel gating. In Section II we introduce the 
model, in Section III we consider a two-state (two-well) 
model which may represent an ion channel with two sta- 
ble conformational states only: one open (where the ions 
can go through the cell membrane) and one closed. In 
Section IV we analyze a more realistic model of gat- 
ing kinetics for a Shaker potassium channel developed 
by Bezanilla et al. [17] based on patch-clamping exper- 
iments with ionic and gating currents. This model was 
also used in [T3] to show the existence of the focusing 
effect. 



II. THE MODEL 

In the model we use, the potential minima are the pos- 
sible states of the ion channel, while the barriers heights 
between two neighboring states are chosen so to repro- 
duce the correct transition rates. It is is described by the 
overdamped Langevin equation 



7 i(i) = -u'(x) + y/2jk B T£(t) 



(1) 



where £(t) is a white noise with (£(t)) = and 
(£(i)£(0)) = S(t), u(x) is the potential landscape, 7 is 
the friction coefficient and T is the temperature of the 
bath. 

We will consider the effect of a time dependent stochas- 
tic external potential on the out of equilibrium stationary 
distributions of our model. When an external voltage, 
v e xt(t), is applied across the cell membrane, the particle 
behaves as if it had a charge valence, z, and the Langevin 
equation becomes (see [H]), 

7±(t) = -u'(x) + zv ext {t)/X+y/2 7 k B T^t), (2) 

where the length scale, A, is the cell membrane thickness. 
Using the rescaled variables 



X = 



A 



T= ^A^' 



T/ U\ V ext(t) TT , s U{x) 

V ext {t) = , U{x) = —- , 

k B l k B l 



(3) 



we finally get the Langevin equation in the dimensionless 
form 



X(t) = -U'(X) + V ext (r) + V2£(r). 



(4) 



In this equation we see that the membrane potential, 
V ex t, changes the potential shape in a linear way: 



U{X) -► U(X) - V ext X. 



(5) 



The specific choice of the external stochastic potential 
is a a dichotomous noise : 



V ext (T) = V + V DN (T), 



(6) 



where Vq is a constant term and Vdn{t) represents a 
is time-dependent random and asymmetric switching be- 
tween to fixed values of the potential, V±,: 



V + = 



D fl + e 



T v \ 1 - 6 



V- = -< 



D fl-e 



t„ V 1 + e 



(7) 



with transition probabilities, w + (from plus to minus 
state), and u>_ (from minus to plus state), 

w+ = (1 + e)/2r v , w- = (l- e)/2r„. (8) 

It is easy to show that such DN has a null time average 

(V dn (t))=0, (9) 

and correlation function 

(V dn (t)V dn {0)) = {D/ Tv )exp(-T/T v ), 

so that it can be fully characterized by three dimension- 
less parameters: correlation time t v , intensity or ampli- 
tude D/t v and asymmetry e (with — 1 < e < 1). 



III. FOCUSING IN A DOUBLE- WELL 
POTENTIAL 

In this Section we consider a two state model for an 
ion channel, represeting open/close configurations. Even 
if it is too simple to adequately represent the function 
of real channels, It is commonly used as toy models, to 
demonstrate various phenomena (T^l riEH2Tj . A two-state 
model can also be understood as a sub-model of a larger 
system, see Section IV. 

The double well potential landscape here considered is: 



u(x) 



1 



1 



-bx 4 



(10) 



which has two minima in ±i m = ±v/ a/b, and a central 
barrier with height Au = a 2 / (46); see Fig. [l] 

Due to thermal fluctuations, the particle is forced 
to jump between the two wells with a rate that for 
Au/k B T 3> 1, is given by the Kramers formula [T]: 



w K 



2-7T7 



exp 



Au 
k^T 



(11) 



where Uq = \u"(x m )\, ui 2 — \u"(xb)\ (evaluated at the 
top Xb of the central barrier). 

In this Section we consider the effect of a DN pertur- 
bation with Vo = 0, see Eq.(j6|. Depending on the state 
of the DN, there are two possible configurations for the 
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potential: a (— ) configuration and a (+) configuration, 
as depicted in Fig. [TJ For each configuration there are 
two different thermal transition rates: for example, in 
the (— ) configuration w L (— ) stands for the transitions 
from the left to the right well and w R {— ) for the transi- 
tions from the right to the left well; and similarly in the 
(+) configuration. 

In order to study the stationary probability distribu- 
tion one has to solve the stochastic equation, Eq.(|4|. In 
our case the Langevin equation Q has been solved nu- 
merically with an order 1.5 strong Ito-Taylor scheme (in 
explicit form), described by Platen and Wagner in [22] . 
Starting with different initial conditions (random distri- 
bution of initial position) and integrating in time the 
Langevin equation, it is possible to compute the proba- 
bility Pl/r{ t ) to be found at the time r in each potential 
well. After a transient time, we verified that this proba- 
bility reaches a stationary value, P^ , dependent on the 
DN parameters. 

Because of the DN, the potential switches between two 
configurations, so it is possible to define four thermal 
transition times in the system: t l ' r (+/—), where, for 
example, t l (+) is the mean time after which a particle 
jumps from the left to the right well when the potential 
is in the (+) configuration (that is when Vdn = V+). 
Note that in this Section we consider a regime in which 
both the DN amplitude and the thermal fluctuations are 
not too large, compared to the height of the potential 
barrier. Small thermal fluctuations w.r.t. the potential 
barrier correspond to physiological conditions at room 
temperature for ion channels, [13]. Moreover also mod- 
erate DN amplitude is a realistic regime since a large 
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FIG. 1. (Color online) The unperturbed double-well potential 
u{x) of Eq. ( |10[ | and the two configurations of the potential 
when it is modified by the dichotomous noise, as indicated 
in the legend. Potential parameters are a = 10 and 6=1, 
while the DN shown refers to: D/t v = 0.71, k B T = 0.15Au, 
e = 0.8. 




FIG. 2. (Color online) Equilibrium probability in the left 
well versus the DN correlation time r v . Full (black) cir- 
cles have been obtained from the Langevin equation, while 
full (red) line has been obtained from a master equation ap- 
proach. Thermal transition times t l ' r (±) are shown as ar- 
rows (except from r R (+), which is found to be of order 10 6 
and is out of scale). The horizontal upper dashed line rep- 
resents the limit p™ ean ^ computed with Eq. ( [Tit, while the 
dot-dashed horizontal lower line isP/ oc , see Eq7J27j) . Here we 
have e = 0.8, k B T = 0.15Ait, D/r v = 0.71, as in Fig.[TJ 



amplitude would destroy the ion-channcls. Under these 
assumptions, the thermal transition times can be com- 



puted through the Kramers rates, w ' (±) see Eq.(ll|, 



for the perturbed potential, so that we have: 

t l ' r {±) = l/w L/R {±). (12) 

We should also take into account that there are two more 
characteristic times, given by the DN, namely, 

±_ J__ _A_ 
Tv ~ w ± ~ (l±ef v ' 



(13) 



where w± are given by Eq.(|8|. Note that r^f are an 
estimate of the times spent in each configuration between 
two DN switchings. 

The results of Langevin equation are shown in Fig. [2] 
at fixed DN intensity, D/t v , asymmetry, e, and varying 
the DN correlation time, t v . Each data in the graph rep- 
resents the left well equilibrium probability, P^ q . As one 
can see there is a minimum of P^ q , which means that 
the 90% of probability is concentrate in the right well, 
at some peculiar value of t v . This is the nonequilibrium 
kinetic focusing, which should not be confused with the 
trivial focusing which occurs at large t v , where P^ q ~ 0.8. 
Indeed this focusing for large t v can be understood based 
on equilibrium consideration. Indeed, when the DN tran- 
sition times, t~, are much larger than the thermal transi- 
tion times, t L ' R , we can assume that the system reaches 
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separately the thermal equilibrium in each potential con- 
figuration. 

Therefore if P L / R (±) are the thermal Gibbs probabil- 
ities for each well in the (+) and (— ) configuration, then 
it is reasonable to assume that, for — > oo, the proba- 



bility Pf 9 approaches the mean value 



(Tv + T V ) 



mi a ) i 

L 



TDmean 
L 



(14) 

is shown as the dashed horizontal line in Fig. [2] 
and agrees very well with numerical data. 

On the other side, in the limit of fast switching of the 
DN (t v — > 0), the probability distribution becomes sym- 
metric: 



Pl\P%^ 1/2, 



(15) 



which is the value obtained with our numerical simula- 
tions of the Langevin equation, see Fig. [2] The limit 
t v — > can be understood using the Langevin equation, 
by assuming, for simplicity, a periodic, rather than a 
stochastic switching, between the two potential config- 
urations. Indeed, after the first step, assuming we are in 
the (+) configuration, we can write, see Eq.Q: 

(AX)+ = [-U'(X) + V + + V2Z]t+, (16) 

and after the second step 

(AX)- = [-U'(X) + V- + V2C]r-, (17) 

neglecting the terms proportional to t+t~ . The total 
variation after a time t+ + r,7 is then: 



AX = (AX)+ + (AX)" 

= \-U\x) + V2£\(t+ + T~) + (v + r+ + «_t-). 



(18) 



and, dividing by + t v , we get the Langevin equation 



for a fast oscillating potential: 



X = -U'{X) + V2£- 



y- 



(19) 



The last term corresponds to the time average of the DN, 
which is zero (see Eq. Q): this proves that for fast oscil- 
lations of the perturbing potential the Langevin equation 
does not feel the perturbation, namely it behaves sym- 
metrically, as anticipated. 

In the intermediate regime we have a non trivial result: 
from Fig. [2J it is clear that the left well probability P^ q 
has a minimum (i.e. there is a "focusing" in the right 
well) for a particular choice of the DN correlation time 
t v . The focusing effect is an interesting result, since the 
applied DN has zero mean (Eq.([9]), and nevertheless it 
can induce stationary distributions which focus the sys- 
tem in one well. 



In order to understand what are the conditions for the 
focusing, we can attempt a simple explanation through 
a master equation approach: 



P = W(r)P, 



(20) 



where P(r) = (Pl{t),Pr(t)) and the time dependent 
transition matrix W(t), which switches between the two 
configurations : 



W(±) 



-w L {±) w R (±) 
w L {±) -w R {±) 



(21) 



with the same rates as in Eq. In Fig. [2] we compare 
the results obtained through the Langevin equation with 
those obtained via the master equation, Eq.(|2"0|). As one 
can see, there is a good agreement between the two ap- 
proaches when the DN correlation time r v is larger than 
the value at which the best focusing occurs, while, when 
t v — Y the master equation gives a complete different 
result. 

This is a limitation of the master equation formalism, 
for which we can give an explanation by assuming, as 
was done above for the Langevin equation, a periodic 
switching. At the initial time the matrix is in the W(+) 
configuration, and after a time rj" it changes into the 
W{— ) configuration, where it stays for a time t~; and so 
on. Thus assuming that r+ and t~ are both infinitesimal, 
after the first step we can write: 



where (AP) + 
step we have: 



(AP)+ = W(+)t+P(0), (22) 
= P(t+)— P(0). Similarly, after the second 



(AP)" = W(-) t- P(r+) = W(-) T~ P(0), (23) 

where the last term holds neglecting the terms propor- 
tional to t^t~ . The total variation after a time r+ +t~ 
is: 



AP = (AP)+ + (AP)~ 

= [W(+)t+ + W(-)t-]P(0), 



(24) 



and, dividing both sides by t+ + t v , we get a system of 
differential equations with an averaged transition matrix: 



(W) 



W(+)r+ + W(-)t- 



(25) 



•v T i v 



This averaged matrix is different from the unperturbed 
transition matrix, that we would expect for t v — > 0, and 
determines the following equation for the probabilities: 



Pl = (W)hPl + (W) 12 Pr 

= ((W) 11 ~(W) 12 )P L + (W} 12 . 

and the stationary solutions P^ — 0, is given by: 



P'r 



(W) 12 -{W)i 



p. 



foe 



(26) 



(27) 
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sition time, t l (+) for different asymmetry e, the inten- 
sity D/t v and the temperature T. Data are shown in 
Fig. [3j Each value of r+ min has been computed from 
a graph like Fig. [2j by making a quadratic fit of the 
function P^ q (log(T v )) around its minimum. In the same 
figure each curve at fixed D/t v and fcsT, has been ob- 
tained varying the asymmetry parameter e. We have the 
remarkable result that r+ min depends on the amplitude 
of noise and temperature only through r L (+). Interest- 
ingly, the relationship between T+ min and t l (+) is more 
complicated than what is suggested in Eq. ( 28 ) . Indeed in 



Fig.[3]the dashed line indicates a square root dependence, 
so that we have: 



v mm 



(29) 



FIG. 3. (Color online) Correlation time r„ min , for which we 
observe focusing, vs the transition time t l (+). Each data set 
refers to simulations with different e but same intensity D jr v 
and temperature T as indicated in the legend. The dashed 
line indicates a square root dependence. Error bars describe 
the variation of r+ min for which changes by 0.1% from 
its minimum value. 



Therefore this average transition matrix in general does 
not have a symmetric solution. In Fig. [2] the value of 
P[ oc is shown to agree with the master equation solution 
in the limit t v — > 0, and, interestingly, it is close to the 
focusing value of P^ q at the minimum. This coincidence 
will be used below to determine the dependence of the 
focusing intensity on the parameters of the model. 

An intuitive explanation of the failure of the master 
equation approach for t v —¥ 0, is that the system does not 
have time to equilibrate in each well before the switching 
and thus the two-state model becomes inappropriate. 

The fact that a master equation approach cannot ex- 
plain the focusing effect which occurs for fast switching 
of the DN, shows that the focusing is highly non trivial, 
being a truly non-equilibrium effect. For a better under- 
stand this focusing we need to consider the interplay of 
several time scales. Physically, we can expect focusing 
when r+ (mean time for which the potential is in (+) 
configuration) is of the same order than t l (+) (i.e. the 
time the particles need to jump towards the right well in 
the (+) configuration) and, at the same time, r~ should 
be much smaller than the time t r (— ) the particles need 
to jump towards the left well. Roughly speaking we can 
thus write the focusing conditions at the right well as, 



r-«r«(-)«r i (-). 



(28) 



We have no theoretical arguments suggesting the scaling 
law above, and more investigation is required. 

Up to now, we have only studied the probability dis- 
tribution induced by the DN as a function of t v . Now we 
intend to analyze the dependence of the focusing on the 
other parameters, such as e and ksT . In Fig. [4] (upper 
panel) we show the left well probability at the focusing, 
P^ q vs the DN asymmetry parameter, e. We observed 
that the focusing decreases as the asymmetry decreases 
(e — > 0). This can be explained by notecing that the dis- 
tance between the two thermal transition times t l (+) 
and t r (— ) decreases as the asymmetry decreases, there- 
fore the second focusing condition in Eq.(28) is not met. 
In Fig. [3] (lower panel) we show the left well probability 
at the focusing, P^ q vs the temperature. Obviously the 
temperature destroys the focusing effect. 

From Fig. [2] we noticed that the master equations fail 
to fit the data when t v is smaller than the value at which 
a minimum is found. After that point, the equilibrium 
probability computed with the maste r eq uations reaches 
the limiting value, P[ or , given in Eq.(27), which is close 
to P?_i evaluated at the focusing. Thus we can use 



Eq.(27) to analytically estimate the value of the focus- 
ing as a function of the parameters. The result is shown 
in both panels of Fig. [4] as dashed curves, which are in 
good agreement with our numerical results, as long as 
the asymmetry, e, and the temperature, kgT, are not 
too large. Indeed in these cases the Kramers formula 
( 1 1 1 is no longer applicable since the energy barriers are 



not well defined. 

The results of our investigations on a double well po- 
tential can be summarized as follows: 

• a stochastic perturbing potential with zero time av- 
erage can induce an unbalanced stationary proba- 
bility distribution between the two wells. 

• conditions for the right-well foc usin g are: r+ ~ 



We now study how the switching time, T+ min , at which 
a minimum of Pf q is found, depends on the thermal tran- 



t l (+) and r~ < t ( — ) , see Eq.([28|). Similarly for 
the left-well focusing. 

• At fixed DN intensity the focusing can be improved 
by increasing the asymmetry e or decreasing the 
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FIG. 4. (Color online) Upper panel: P^ q at focusing (black 
dots) is plotted vs e for D/r v = 0.71, ksT = 0.15Am. Lower 
panel: P^ q at focusing (black dots) is plotted vs UbT/ Au for 
D/t v = 10(/cbT) 2 , e = 0.8. In both panels the dashed lines 
represent the theoretical predictions obtained from the master 
equation, see text. 



temperature. Analytical estimation of the focusing 
intensity is given in Eq.(27|. 



• The master equation approach completely fails for 
fast DN perturbations (small r„) and is unable to 
describe the nonequilibrium kinetic focusing. 



IV. EIGHT- WELL POTENTIAL: A MORE 
REALISTIC MODEL FOR ION CHANNEL 
GATING 

In the previous Section we showed that it is possi- 
ble to modify the stationary probability distribution in 
a double-well potential, and particularly to increase the 
probability to find the system in one well, by introduc- 
ing a stochastic perturbation. Although a potential with 
only two states can be a useful starting point, a multi- 
well potential is a more realistic model for ion channels. 
As an example, we consider, as in [T3], the Shaker K + 
channel, for which there exists a simple kinetic model, 
proposed by Bezanilla, Perozo and Stefani (BPS) [T7j: in 
this model, there are eight states, or conformations, and 
the system can jump between them. All the states from 
1 to 7 correspond to closed conformations, and the state 
8 is the only one for which the ion channel is open (which 
means that the ions can go through the channel) . The al- 
lowed transition rates for the Shaker K + ion channel are 
known from experimental measurements. Since transi- 



tions occur only between nearest neighboring states, this 
can be represented by an eight-well potential landscape 
u(x), as depicted in Fig. [H] The eight potential minima 
are the possible states of the channel, and the barriers 
heights are chosen so to reproduce the correct exponen- 
tial transition rates. 




FIG. 5. The potential u used by Millonas and Chialvo in [14] 
as a model of the Shaker K + channel. For the definition of 
the rescaled potential U(X), see Eq. |3|. In this Figure we 
assumed T = 3007^. 

In the following we present some preliminary discussion 
on the effect of stochastic noise on a multiwell potential 
landscape, leaving a deeper analysis for future work. 

The problem of focusing in this eight-well model has 
been already considered by Millonas and Chialvo [14j . 
Based on a probabilistic treatment of the Langevin equa- 
tion Q, Millonas and Chialvo claimed that, for very low 
temperature and large amplitude of the DN, it is possible 
to induce stationary probability distributions located in 
a chosen potential well: they called this effect nonequilib- 
rium kinetic focusing. These two conditions (low temper- 
ature and large DN) are not physiological, since the ion 
channel would be destroyed under these conditions. For 
this reason, we will first investigate whether a nonequilib- 
rium kinetic focusing can occur in an eight-well potential 
in physiological conditions (room temperature and mod- 
erate DN amplitude) as it was considered in the double 
well case in the previous Section. We will also compare 
the analytical predictions of Ref. jUj for low temperature 
and large DN intensity with our numerical simulations. 



A. Room temperature and low DN intensity 

Our starting point is the Langevin equation Q, and 
the analysis of the results will concern the equilibrium 
probabilities P^ q in each well, denned as in the previous 
Sections. 
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In Fig.[6]we report our simulations concerning the equi- 
librium probability distribution P^ q as a function of t v 
for each well, at fixed asymmetry e, intensity D/t v and 
mean external potential Vo- Here we can see that some 
kind of focusing occurs only close to the potential edges 
(states 1, 2 or 7, 8), while the equilibrium probabilities of 
the central wells are almost independent of t v . 




FIG. 6. (Color online) Density plot of the equilibrium proba- 
bilities P^ q in the i-th well versus the DN correlation time r v . 
Darker color intensities stand for higher probabilities. Data 
are: D/t v = 1, e = 0.8, Vo = {V ext } = (-45 mV)/fc s T, and 
T = 300 K. 

From the same figure it is also possible to recover the 
usual observations made in the double-well case about 
the limits of very small and very large r v . We checked 
that when t v — »• 0, the DN averages to zero, so that the 
probability distribution the thermal Gibbs distribution 
for V ext = Vo- On the other hand, for t v ^> 1 the proba- 
bilities can be evaluated, as in Eq.(14), by 
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where the probabilities -Pi(±) are the thermal Gibbs 
probabilities for each well in the (+) and (— ) configu- 
ration. 

The non trivial focusing in the intermediate r„ regime, 
resembles the one found in the double well case. Indeed 



the conditions given in Eq. ( 28 ) seem to work also in the 



multi-well case if we consider the edge wells 1 , 2 as a left 
well; the edge wells 7, 8 as the right well and the central 
wells as an effective barrier between them. Nevertheless 
a more detailed analysis is necessary to understand how 
to relate the focusing in the double well with that in the 
multi-well case. 

We also considered many other situations, with differ- 
ent values of the DN asymmetry parameter e and of the 
DN intensity D/t v : in all cases, we could only observe 
some probability increase in the wells on the edges, while, 
to the best of our effort, we were not able to obtain focus- 
ing in the central wells. Such a result would be desirable 
since it would allow a complete control of the ion channel 
dynamics under a stochastic perturbations. 



B. Low temperature and large DN intensity 

An analytical solution (Eq.(7) in Ref. [14]) for the sta- 
tionary probability distribution induced by a DN, were 
derived under the following assumptions: 

• the amplitude of the driving, must be larger than 
the thermal fluctuations, 

• and larger than the maximal potential jump, 
y/fifc>8up\U'\; 

• the temperature should be very small (that is T — > 
0). 

In the above Section, we did not take these conditions 
into account, because we considered them not physio- 
logically feasible, nevertheless in this Section we take 
this regime into consideration in order to understand 
the range of validity of the analytical prediction given 
in Ref. [H]). 

In Fig. [7] we compare our numerical results (hys- 
togram) with the theoretical results (dots) given in 
Ref. [14] for different temperatures while fixing all the 
other parameters. As we can see, while at low temper- 
ature (panel a) and b)) there is a good agreement, on 
increasing the temperature (panel c) and d)) the the- 
oretical distribution starts to fail even for temperature 
well below room temperature (T — 40 K). Interestingly 
we notice a good focusing in the central wells: in panel 
(c) Langevin equation gives approximately 90 % of prob- 
ability in the wells 5 and 6), even when the theoretical 
equation of Ref. [U] does not work. The presence of some 
focusing in this case seems to indicate that the conditions 
given in Ref. [14] , could be too restrictive in order to have 
focusing. 



V. CONCLUSIONS 

In this work we considered two different kinetic mod- 
els of ion channels subjected to dichotomous noise per- 
turbation, with zero time average. The first is a sim- 
ple model consisting of two wells separated by an energy 
barrier, while the second is the physiologically relevant 
eight-well model proposed by Bezanilla et al. [17]. In 
both cases an highly non trivial non equilibrium focus- 
ing has been found under physiological conditions (room 
temperature and moderate external perturbation inten- 
sity), which cannot be described by a master equation 
approach, but the analysis of the full Langevin equation 
turns out to be necessary. Conditions for the focusing 
and dependence on the external perturbation parameters 
have been given for the double well case. In the eight- 
well case more analysis is needed in order to address this 
question. 

We also analyzed the results of Millonas and Chialvo 
in |14j , showing that the nonequilibrium kinetic focusing 
does not survive in physiological conditions. 
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In perspective our analysis leaves still open the ques- 
tion of whether it is possible to focus an ion channel in 
any well, thus completely controlling its dynamics, un- 
der physiological conditions. Nevertheless we have shown 
both that focusing is possible in physiological conditions 
and also we gave a preliminar evidence that focusing in 
the central wells of a multi-well potential is also possible. 



FIG. 7. (Color online) Probability distributions at dif- 
ferent temperature computed with the Langevin equation 
(green bars) and theoretically (black diamonds) with Eq.(7) 
in Ref. [T5]. The parameters are: D/t v = 10 5 , e = 0.6, 
t v = 10" 4 , v = -142 mV/k B T. Temperature is : T = 20K 
(a), T = 25K (b), T = 30K (c), T = 40K (d). 



[1] L. Gammaitoni, P. Hanggi, P. Jung and F. Marchesoni, 
Stochastic Resonance, Rev. Mod. Phys. Vol. 70, No. 1 
(1998) 

[2] J. Kallunki, M. Dube, T. Ala-Nissila, J. Phys:Cond. 

Matt. 11, 9841 (1999) 
[3] S.M. Bezrukov, I. Vodyanoy, Nature 378, 362 (1995) 
[4] P. Hanggi, Chem. Phys. Chem. 3, 285 (2002) 
[5] R.K. Adair, PNAS 100, 12099 (2003) 
[6] I. Goychuk, P. Hanggi, J.L. Vega, S. Miret-Artes, Phys. 

Rev. E 71, 061906 (2005) 
[7] I. Goychuk, P. Hanggi, Phys. Rev. Lett. 91, 070601 

(2003) 

[8] K. Wiesenfeld, F. Jaramillo, Chaos 8, 539 (2001) 
[9] M.J. Chacron, G.W. Slater. Phys. Rev. E 56, 3446 (1997) 
[10] C.R. Doering, J.C. Gadoua, Phys. Rev. Lett. 69, 2318 
(1992) 

[11] O. Flomenbom, J. Klafter, Phys. Rev. E 69, 051109 
(2004) 

[12] G.J. Schmid, P. Reimann, P. Hanggi, J. Chem. Phys. 



Ill, 3349 (1999) 
[13] K. Lee, W. Sung, Phys. Rev. E 60, 4681 (1999) 
[14] M. M. Millonas, D. R. Chialvo, Phys. Rev. Lett. 76, 550 

(1996) 

[15] A. Kargol, K. Kabza, Phys. Biol. 5 (2008) 026003 

[16] B. Hille, Ionic channels of excitable membranes, Sinauer 

Associates Inc. (1992) 
[17] F. Bezanilla, E. Perozo, E. Stefani, Biophys. J. 66, 1011 

(1994) 

[18] M.C. Mahato, A.M. Jayannavar, Physica A 248, 138 
(1996) 

[19] P.K. Ghosh, D.S. Ray, J. Chem. Phys. 125, 124102 (2006) 
[20] E. DiCera, J. Chem. Phys. 95, 5082 (1991) 
[21] M.C. Mahato, A.M. Jayannavar, Phys. Lett. A 209, 21 
(1995) 

[22] P. E. Kloeden, E. Platen, Numerical Solution of Stochas- 
tic Differential Equations, Springer (1995) 



